arXiv:1508.00567vl [astro-ph.HE] 3 Aug 2015 


Mon. Not. R. Astron. Soc. 000, 1-11 (2011) 


Printed 5 August 2015 


(MN [AT E X style file v2.2) 


New constraints on the structure and dynamics of black hole jets 

William J. Potter* and Garret Cotter 

Oxford Astrophysics. Denys Wilkinson Building, Keble Road, Oxford, OX1 3RH, United Kingdom 


ABSTRACT 

Accreting black holes produce powerful relativistic plasma jets which emit radiation across 
all observable wavelengths but the details of the initial acceleration and confinement of the jet 
are uncertain. We apply an innovative new model that allows us to determine key properties of 
the acceleration zone via multi-frequency observations. The central component of the model 
is a relativistic steady-state fluid flow, and the emission from physically distinct regions can be 
seen to contribute to different energy bands in the overall spectrum. By fitting with unprece¬ 
dented accuracy to 42 simultaneous multiwavelength blazar spectra we are able to constrain 
the location of the brightest synchrotron emitting region, and show that there must be a linear 
relation between the jet power and the radius of the brightest region of the jet. We also find 
a correlation between the length of the accelerating region and the maximum bulk Lorentz 
factor of the jet and find evidence for a bimodal distribution of accretion rates in the blazar 
population. This allows us to put constraints on the basic dynamical and structural properties 
of blazar jets and to understand the underlying physical differences which result in the blazar 
sequence. 

Key words: Galaxies: jets, galaxies: active, radiation mechanisms: non-thermal, radio contin¬ 
uum: galaxies, gamma-rays: galaxies. 


1 INTRODUCTION 

Accreting black holes can launch powerful relativistic plasma jets 
which are thought to be produced by the interaction of a rotating 
black hole with the magnetic field coupled to the accreting plasma. 
The rotation of the black hole induces rotation of inertial frames sur¬ 
rounding the black hole which twists the magnetic field close to the 
event horizon. This process transfers energy from the black hole’s 
rotation into the energy contained in the magnetic field (Blandford 
& Znajek 1977). The enhanced pressure of the magnetic field is 
strong enough to drive and accelerate a plasma jet, however the de¬ 
tails of the acceleration and confinement of the jet are uncertain 
(Begelman et al. 1984, Vlahakis & Konigl 2003, Hawley & Krolik 
2006 Beskin & Nokhrina 2006, McKinney 2006, Komissarov et al. 
2007, Tchekhovskoy et al. 2010 and McKinney et al. 2012). Radio 
observations of the jet in M87 show that it starts with a parabolic 
shape, transitioning to a conical jet at 10 5 Schwarzschild radii, r., 
(Asada & Nakamura 2012). These observations are consistent with 
general relativistic magnetohydrodynamic (GRMHD) simulations 
in which the jet starts with a magnetically dominated parabolic ac¬ 
celerating base and transitions to a ballistic conical jet. This tran¬ 
sition occurs when the jet plasma comes into approximate equipar- 
tition between the energy contained in the magnetic field and par¬ 
ticles. Here we use the first jet emission model which takes into 
account these state-of-the-art results from observation and theory 
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in order to constrain the basic structure and dynamics of jets by 
fitting to a large sample of blazar spectra. 

The most luminous class of supermassive black hole jets are 
blazars: jets which are pointing towards Earth so their emission 
is strongly Doppler-boosted. Evidence for their Doppler-boosting 
comes front the observed superluminal motion within their jets and 
short-timescale variability (Urry & Padovani 1995, Ulrich et al. 
1997, Jorstad et al. 2001 and Abdo et al. 2010). Non-thermal elec¬ 
trons are accelerated within the jets by processes which are cur¬ 
rently not well understood (most likely a combination of shocks 
and magnetic reconnection e.g. Summerlin & Baring 2012; Cerutti 
et al. 2012; Nishikawa et al. 2013) and these electrons emit radia¬ 
tion across all observable wavelengths via synchrotron and inverse- 
Compton mechanisms. Blazars are particularly important to our 
understanding of jet properties because their emission is strongly 
Doppler-boosted and so the entire multiwavelength jet spectrum is 
clearly observed above the galactic emission. This makes it possible 
to constrain basic jet properties by using realistic jet emission mod¬ 
els to fit to observed spectra. Previous models of blazar jet emission 
have focused either on understanding the high energy gamma-ray 
emission by modelling the jet as one or two fixed radius spherical 
blobs (e.g. Kirk et al. 1998, Li & Kusunose 2000 Bottcher & Chi- 
ang 2002. Tsang & Kirk 2007, Dermer et al. 2009, Tavecchio et al. 
2011, Bottcher et al. 2013, The MAGIC Collaboration et al. 2014, 
Petropoulou 2014 and Pacciani et al. 2014), or using extended coni¬ 
cal models to fit primarily to the flat radio spectrum (e.g. Blandford 
& Konigl 1979, Mufson et al. 1984, Ghisellini et al. 1985, Kaiser 
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et al. 2000, Markoff et al. 2001, Spada et al. 2001, Pe’er & Casella 
2009, Jamil et al. 2010 and Vila et al. 2011). 

In this paper we use one of the most sophisticated models for 
blazar jet emission currently available (Potter & Cotter 2012, Potter 
& Cotter 2013a, Potter & Cotter 2013b and Potter & Cotter 2013c) 
to fit to the entire sample of Fermi blazars from Abdo et al. (2010) 
with simultaneous multiwavelength observations and redshifts. Be¬ 
cause this model takes into account the extended parabolic to con¬ 
ical geometry and acceleration of the jet fluid, we show that we 
are able to place meaningful constraints on the structure and dy¬ 
namics of supermassive black hole jets by fitting to the blazar spec¬ 
tra across all observable wavelengths, with unprecedented accuracy. 
The structure of the paper is as follows: we start by introducing and 
explaining the assumptions of our jet model, we then show the re¬ 
sults of fitting the model to the observed spectra and the constraints 
we obtain. Finally, we discuss these results in the context of cur¬ 
rent ideas on magnetic acceleration, the blazar sequence and AGN 
unification. 


2 JET MODEL 

Our model is motivated by the recent results from observations and 
simulations. The jet is modelled by a ID time-independent relativis¬ 
tic fluid flow with a variable shape and bulk Lorentz factor (Potter 
& Cotter 2013a). The total relativistic energy of the plasma is con¬ 
served via the equation for energy-momentum 
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where is the total energy-momentum tensor of the jet plasma 
which can be decomposed into magnetic and particle energy den¬ 
sities, and also a cumulative energy loss term which we include 
to conserve the total energy along the jet. The components of the 
energy-momentum tensor in the fluid rest frame, indicated by a 
prime, are given by 
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where we have set c = 1, x is the distance along the jet axis in 
the lab frame, B' is the rest frame magnetic field strength, E e is the 
electron energy, n e the electron energy distribution and T > s , y nc h( a: ) + 
P{c( x ) + -Pad the sum of the synchrotron, inverse-Compton and 
adiabatic losses per unit length in the fluid rest frame (for a detailed 
calculation of these loss terms see sections 3-6 in Potter & Cotter 
2012 and 5-6 in Potter & Cotter 2013a). Making the assumption 
that the plasma is locally homogeneous and isotropic perpendicular 
to the jet axis, the energy-momentum tensor simplifies to that of a 
relativistic perfect fluid, P = p/3, in the fluid rest frame. 
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from the rest frame to the lab frame (and assuming /? = v/c « 1 ) 
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where 7 buik is the bulk Lorentz factor of the jet plasma. The two 
factors of the Lorentz factor in the relation of the lab frame energy 
density to the rest frame energy density, p oc p’"f 2 , can be intu¬ 
itively understood as one being due to an increase in energy by the 
Lorentz boost and one due to the increase in density because of the 
length contraction in the x-direction. Integrating (1) and using the 
4-dimensional divergence theorem (see for example Page & Thorne 
1974) we find 
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where, cPV = \Z\g\dtdxdRd(/>, is the invariant 4-volume, g, is the 
determinant of the metric tensor and \f\g\ = R, using cylindrical 
coordinates in Minkowski space (appropriate for the strongly radi¬ 
ating sections of the jet which occur many Schwarzschild radii from 
the black hole). The last three 3-surface integrals containing an inte¬ 
gral over time are equal to zero due to the time-independence of the 
model and we have used the assumed radial and azimuthal symme¬ 
try of the jet plasma. The only non-zero components of (5) occur for 
v — 0 or 1 , in both cases this gives us our equation for conservation 
of energy-momentum 

^7buik(a:) 2 7ri? 2 (a:)p , (x)^ = 0. (6) 

The particle flux J M is conserved along the jet by the equation 

V M x) = 0, r{E e ) = n e (E e ,x)U >1 (x). (7) 


where U^/x) = 7 (x)(l, (3(x), 0,0) is the jet fluid 4-velocity, 
P{x) = v{x)/c , v(x ) is the jet speed and n' e the electron number 
density in the rest frame. Integrating (7) and using the divergence 
theorem as before we find 

j V M J Al dV = J = ^{■KR 2 {x)n' e (x)U°(x)) = 0 , 

( 8 ) 

where again the three 3-surface integrals which contain an integral 
over time vanish due to the time-independence of the model. These 
equations ensure that the total energy is conserved along the jet and 
naturally take into account the magnetic energy required to accel¬ 
erate the jet and the internal energy gained when the jet decelerates 
via interactions with its environment. To illustrate this we shall ex¬ 
plicitly show that in the case of an accelerating jet with minimal 
radiative and adiabatic losses the magnetic energy is converted into 
bulk kinetic energy of the plasma. Equations 6 and 8 have solutions 
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where p' = T ,ol) is the total rest frame energy density. Using 
a Lorentz transformation to convert the energy momentum tensor 
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Transition region. Jet transitions from parabolic to conical. 
Plasma first comes into equipartition and magnetic 
acceleration ceases to be efficient. Dominates optically thin 
synchrotron and SSC emission. 



Figure 1 . Jet schematic. 


The electron energy density is given by (2) and we decompose p' 
into the rest frame magnetic and particle energy densities. We shall 
temporarily neglect the energy loss term in this calculation for clar¬ 
ity. 
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Since the total rest frame energy density is decreasing as 7 j^ k , 
whereas, the rest frame particle energy density is only decreasing as 
7buik’ w h en the jet accelerates and the bulk Lorentz factor increases 
the jet plasma will become more and more particle dominated. This 
is because magnetic energy is being converted into particle energy 
as the bulk jet plasma is accelerated by magnetic forces and even¬ 
tually, when the magnetic energy becomes smaller than the bulk 
particle energy, this acceleration ceases to be efficient. This demon¬ 
strates the benefit of using the relativistic fluid equations 6 and 8 as 
the basis of the jet model. 


2.1 Model assumptions 

The best current constraints on the shape of a supermassive black 
hole jet come from radio observations of M87 which show a 
parabolic base extending out to a distance, xt = 10 5 r s , where the 
jet transitions to conical (Asada & Nakamura 2012). We use these 
observations as the basis for our jet model choosing our jet shape to 
be 

R(x) = C P ar(r s + x) Apar for x < xt, (13) 

R(x) = C par (r s +XT) ApB ' 1 +(x—X T ) tan(6 l 0 pening) for X > XT, 

(14) 

where r a is the Schwarzschild radius and we estimate A par = 0.58 
and C P ar = 1.49r°' 42 from the observations of M87 (Asada & 
Nakamura 2012), and 9 con is the lab frame conical jet opening an¬ 
gle (we assume that due to relativistic beaming 9 con « l/7buik). 
Hereafter we shall refer to the region where the jet transitions from 
parabolic, accelerating and magnetically dominated to a slowly de¬ 
celerating conical jet in equipartition, as the transition region of the 
jet and we shall use a subscript T to denote quantities measured at 
this transition region. We choose to linearly scale the jet geometry 


from the observations of M87 by using an effective black hole mass, 
M, defined such that the transition from parabolic to conical occurs 
at 10 5 r s (M) = xt for all fits, as observed in M87. This fixes the 
aspect ratio of the transition region radius to transition region dis¬ 
tance ( Rt/xt ~ 1/80) to be fixed to that observed in M87, whilst 
allowing the absolute lengthscales to depend linearly on M. This 
seems a much more reasonable assumption than allowing the jet as¬ 
pect ratio at the transition region to be a function of the absolute 
distance to the transition region (in this case the aspect ratio would 
vary as Rt/xt oc xj Apar using equation 13). The effective black 
hole mass, M, found by fitting the model to a blazar spectrum, rep¬ 
resents the black hole mass that would be required for the transition 
region for this blazar to occur at 10 5 r s as in M87. In section 3.3 we 
discuss observational evidence which indicates that differences in 
the jet geometry are unlikely to be driven primarily by changes in 
the black hole masses and so the reader should consider the effec¬ 
tive black hole mass, M, simply as a convenient parameter to vary 
the distance to, and radius of the transition region. 

In agreement with the observations of M87, GRMHD simu¬ 
lations find that the jet starts with a parabolic magnetically domi¬ 
nated accelerating base (e.g. McKinney 2006). In these simulations 
the jet is accelerated by a magnetic pressure gradient which converts 
magnetic energy into bulk kinetic energy and this process continues 
until the jet plasma approaches equipartition between the magnetic 
field and particle energies (where the magnetic energy is no longer 
large enough to significantly accelerate the plasma). We incorpo¬ 
rate all of this information into our model, a schematic of which 
is shown in Fig. 1. In our model the jet starts with a magnetically 
dominated accelerating parabolic base and we use the relation be¬ 
tween the jet shape and acceleration calculated by Beskin & Nokh- 
rina 2006 and consistent with simulations e.g. McKinney 2006 

7bulk(*) OC X 1/2 . (15) 
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Once efficient acceleration of the jet has ceased and it transitions 
to conical we assume the bulk Lorentz factor slowly decreases due 
to interaction of the jet with its environment i.e. entrainment. We 
choose a logarithmic deceleration so the conversion of bulk kinetic 
energy into in-situ non-thermal particle acceleration is continuous 
along the jet, as suggested by observations (e.g. Jester et al. 2005). 
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In our previous work (Potter & Cotter 2013a) we found that accel¬ 
erating non-thermal electrons in the parabolic base of the jet had 
little observable effect since the bulk Lorentz factor in the base 
is not as large as in the conical jet and so the emission front the 
parabolic regions is dominated by the emission from the faster con¬ 
ical jet, which is more strongly Doppler-boosted. We also found 
that the average equipartition fraction in the base must be very low 
(i.e. almost completely magnetically dominated). This is because 
the timescale for electrons to lose energy by synchrotron emission 
was very short in the high magnetic fields present in the accelerat¬ 
ing region and so the energy required to maintain a significant ob¬ 
servable population of non-thermal electrons was unfeasibly large. 
For these reasons we assume that only a small amount of energy 
is contained in electrons at the base of the jet. We assume that the 
jet comes rapidly into equipartition at the transition region between 
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Figure 2. The energy density of the different external photon sources as a 
function of jet length, calculated for the fit to the FSRQ J1512 in the plasma 
rest frame. The rest frame energy densities account for the Doppler boosting 
and beaming of the external photon sources into the plasma rest frame. The 
acceleration of the bulk Lorentz factor (up to x ~ 10 18 m) of the jet leads 
to the increase in the energy density of the CMB as measured in the plasma 
rest frame and the slow deceleration of the conical jet leads to the decrease 
in the CMB energy density beyond this distance. The synchrotron break 
constrains the transition region of FSRQs to occur outside the expected radii 
of the BLR and dusty torus, at large distances along the jet the CMB is the 
dominant external photon field. We find that inverse-Compton scattering of 
CMB photons at large distances is able to fit very well to the high energy 
X-ray and gamma-ray observations for FSRQs. 


parabolic and conical regions where the jet stops accelerating ef¬ 
ficiently. This rapid particle acceleration at the transition between 
parabolic and conical sections could be caused by a combination of 
a recollimation shock if the jet becomes underpressured relative to 
its environment (e.g. Laing & Bridle 2002, Bromberg & Levinson 
2009) and enhanced magnetic reconnection due to the jet compres¬ 
sion but a detailed consideration of these mechanisms is beyond the 
scope of this paper. In our previous work we found that BL Lac type 
blazars were well modelled by jets which maintain equipartition be¬ 
tween non-thermal electrons and the magnetic field in the conical 
section, whereas flat spectrum radio quasars (FSRQs) favoured an 
equipartition fraction which began in equipartition at the start of 
the conical jet and became progressively more particle dominated 
at large distances. This is consistent with the theoretical interpreta¬ 
tion of jet morphologies by Meier (1999), who suggested that FRII 
jets should be particle dominated at large distances while FRI jets 
remain highly magnetised. 

We have chosen to use a simple electron-positron jet compo¬ 
sition in this work to increase the constraining power of our model. 
We know that high-energy non-thermal electrons exist in jets since 
we see polarised synchrotron emission (Jorstad et al. 2005) and so 
it is sensible to first try to fit observations with an electron-positron 
jet before adding additional components which carry extra free pa¬ 
rameters. Including hadronic emission processes would add addi¬ 
tional free parameters to the model making it more difficult to con¬ 
strain the model. Hadronic emission models currently seem to re¬ 
quire highly super-Eddington jet powers to match to data. For ex¬ 
ample, in Bottcher et al. 2013, the proton kinetic luminosities range 
from 2.1 x 10 39 — 4.4 x 10 43 W, even for a very large 1O 1 O M 0 
black hole, 4.4 x 10 43 W is 350 times the Eddington Luminosity, 
which seems implausibly large. Adding additional components to 


our model such as a cold proton component would simply increase 
the jet power. We can view our model as providing an estimate of 
the power contained in the magnetic and non-thermal electrons in 
the jet and thus an estimate of the lower limit of the jet power. 

We calculate the inverse-Compton emission by Lorentz trans¬ 
forming the external radiation fields into the plasma rest frame and 
using the full Klein-Nishina cross-section. We comprehensively ac¬ 
count for the sources of external radiation as functions of distance 
along the jet (Potter & Cotter 2013a) these include photons from: 
the accretion disc, the broad line region (BLR), the dusty torus, 
the narrow line region (NLR), starlight and the cosmic microwave 
background (CMB). FSRQs are defined observationally as blazars 
in which broadened emission lines from the accretion disc have 
been observed and so for these sources we fit a standard composite 
blackbody thin accretion disc model to the spectra (Shakura & Sun- 
yaev 1973). The distribution of these different external photon fields 
is shown in Figure 2. It is worth emphasising the importance of in¬ 
cluding the CMB photons which we find contribute significantly 
to the gamma-ray emission of FSRQs at large distances. We also 
calculate the amount of 7 -ray absorption due to photon-photon col¬ 
lisions which produce electron positron pairs (see Potter & Cotter 
2013c). 


2.1.1 Acceleration of non-thermal electrons 

The mechanism which accelerates non-thermal electrons in jets is 
not well understood currently, the most popular physical processes 
are acceleration by shocks and magnetic reconnection. These mech¬ 
anisms tend to produce non-thermal electron distributions which 
have a power law component at low energies and a cutoff at high 
energies (Kirk & Schneider 1987; Zenitani & Hoshino 2007; Sum¬ 
merlin & Baring 2012) and we have found that this is compatible 
with observed emission spectra of blazars Potter & Cotter (2013b). 
Our model takes into account the radiative and adiabatic energy 
losses to the electron population travelling along the jet and also 
the acceleration of additional electrons which is required to keep 
the plasma close to equipartition along the jet (as found from ob¬ 
servations e.g. Croston et al. 2005 and Homan et al. 2006). The 
energy for this in-situ acceleration comes from the internal energy 
gained from the deceleration of the jet, where bulk energy of the 
plasma is converted into enhancing the magnetic field strength and 
accelerating non-thermal electrons via shocks, and also from mag¬ 
netic reconnection whereby internal magnetic energy is converted 
into accelerating non-thermal electrons. Since the precise electron 
distributions originating from these processes are far too computa¬ 
tionally expensive to calculate within our model (they require dedi¬ 
cated simulations in their own right) we assume that the initial elec¬ 
tron distribution and the distribution of any additional non-thermal 
electrons accelerated in-situ is given by 

Mniected(z, E e ) = (18) 

where a is the electron energy distribution spectral index. Of 
course, due to the energy losses from the electron population (which 
are a complicated function of the electron energy and jet param¬ 
eters) the total electron distribution will quickly deviate from the 
simple form in (18) as the plasma flows along the jet and radiates. 

This model incorporates current observations and results from 
simulations and theory. We have made the model flexible and rig¬ 
orous but at the same time chosen to make simplifying assumptions 
in order to limit the number of free parameters. The model has only 
12 free parameters which is in fact far fewer than a model with two 
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independent spherical emission zones, however, we argue that our 
model is considerably more realistic and than these models. 


2.1.2 Adiabatic losses in a steady jet 


The external adiabatic energy losses which transfer internal energy 
from a thin cylindrical volume element of jet plasma to the confin¬ 
ing medium occur at a rate 
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where p ext is the pressure of the jet's external environment, E\ nt = 
Ep + Eb is the internal energy contained in the cylindrical vol¬ 
ume element, as measured in the lab frame and Ep and Eb are the 
components of energy contained in particles and magnetic fields. 
We decompose the change in the volume of the cylindrical element, 
dV ( x, t) defined via 


rR(t) 

dV(x,t)= / 2ttR(x,t)dxdR(x,t), (20) 
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into a long-term time average dV(x) and a short-term fluctuat¬ 
ing part dV(x,t). If the time average converges for suitably long 
timescales then this is equivalent to saying that the time average of 
the fluctuations will tend to zero (i.e. at a fixed distance along the jet 
its radius will not systematically expand or shrink with time). Our 
model assumes a jet with a steady long-term average volume and 
so over suitably long timescales the work done by the jet expand¬ 
ing the external medium will equal the work done by the external 
medium compressing the jet. This means that such a steady-state jet 
will not do any net work on the external medium. 

Whilst the lab frame volume occupied by the jet is likely to be 
well-described by a time-independent jet model, clearly individual 
fluid elements will expand in radius as they propagate along the jet. 
In this case the plasma will do work on itself as it expands radially 
by accelerating the bulk radial velocity of the jet (analogous to the 
adiabatic radial expansion and cooling in a supernova explosion). 
The rate of change of relativistic particle energy due to internal adi¬ 
abatic losses is given by (e.g. Kaiser 2006) 

d In Ep (V) 1 

dlnV 3’ ' ; 


Ep(2R) = 0.63Ep(R). (22) 

where V is the volume of a thin cylindrical fluid element and R 
is the jet radius, both measured in the plasma rest frame. For a jet 
which undergoes internal adiabatic losses and has a fixed equipar- 
tition ratio A equ i. The rate of loss of internal energy is modified by 
the store of available magnetic energy to become 

1 Er 
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Where again Ei n t = Eb + Ep. In our model our primary concern 
is to conserve the total energy contained in the jet plasma. It is not 
currently known how the total energy of the plasma is distributed 
between magnetic, particle and bulk kinetic energies along the jet. 
Observations show continuous injection of energy into the accel¬ 
eration of non-thermal electrons along the jet which are observed 
via optical synchrotron emission (Jester et al. 2001 and Jester et al. 


2005). The main candidates for this in-situ acceleration are inter¬ 
nal shocks, which convert bulk kinetic energy into non-thermal par¬ 
ticle energy, and magnetic reconnection which converts magnetic 
energy into particle energy via resistive dissipation. Importantly, 
synchrotron observations favour a continuous acceleration mech¬ 
anism since they do not observe a series of discrete, bright, strong 
shock fronts (Jester et al. 2001). The only processes which lead to 
a net decrease in the total energy contained in the jet plasma are the 
emission of radiation and external adiabatic losses (external adia¬ 
batic losses are zero in a time-independent jet as we have argued 
above). 


Since the energy in the non-thermal electrons requires constant 
replenishment (Jester et al. 2001 and Jester et al. 2005) this leads to 
three possibilities for the source of this energy: resistive dissipa¬ 
tion of the magnetic field, conversion of the bulk kinetic energy of 
both the non-thermal, and possibly an additional component of cold 
particles or the conversion of the radial kinetic energy into internal 
energy. The problem with powering the in-situ particle accelera¬ 
tion with purely magnetic dissipation is that the jet would need to 
be magnetically dominated over the majority of its length and so 
would not remain close to equipartition as suggested by observa¬ 
tions (e.g. Croston et al. 2005). A dominant cold component would 
allow equipartition to be maintained but at the expense of increas¬ 
ing the total jet power to account for this extra source of energy (the 
main effect of this in our current work would be to increase our esti¬ 
mation of the total jet power). The only limitation on this process is 
that it should not require an unreasonably large jet power. The final 
possibility is that the work done by the non-thermal particles ex¬ 
panding the jet (accelerating the bulk radial velocity) is efficiently 
converted back into non-thermal particle energy. This could occur, 
for example, if the jet undergoes recollimation shocks in which an 
initially overpressured jet expands into the external medium until 
it is underpressured, the pressure of the environment then acceler¬ 
ates material radially towards the jet axis where it eventually col¬ 
lides supersonically. This forms a recollimation shock which con¬ 
verts some of the radial kinetic energy back into internal energy via 
non-thermal shock acceleration and magnetic field enhancement. 


In this paper we have chosen this final mechanism in which 
adiabatic losses are replenished by the conversion of radial kinetic 
energy into internal energy. This is because this scenario is the sim¬ 
plest to implement, provides a lower bound estimate of the jet power 
and does not require the introduction of an additional unknown pa¬ 
rameter (the ratio of the energy in cold particles compared to non- 
thermal electrons). We shall investigate the effect of the introduc¬ 
tion of a cold particle component of the fluid flow in a future paper. 


It is worth commenting that the detailed question of whether 
adiabatic losses are really appropriate for a diffuse magnetised and 
non-thermal jet plasma is far from certain. The prolonged existence 
of non-thermal electrons demonstrates that the plasma is collision¬ 
less (particle-particle collisions are infrequent) and so the usual as¬ 
sumptions which go into adiabatic expansion are not fulfilled (the 
system does not move slowly and reversibly between adjacent states 
of thermal equilibrium). Since the number density of non-thermal 
electrons is low, and their energies very large, scattering of the par¬ 
ticles from the field lines is stochastic and is not likely to be well 
described by a fluid approximation. To determine whether adiabatic 
losses apply in their usual form in these extreme conditions might 
be best tackled through particle in cell simulations. 


© 2011 RAS, MNRAS 000, 1-11 



6 William J. Potter and Garret Cotter 


o 

Q. 

c 

o 

~ H—' 

\n 

c 

n3 

i_ 

-*—i 

H— 1 

05 

</) 

=3 


T3 

05 

L. 

H— 1 

CD 

—> 


1 - 


0.1 - 


^ 0.01 - 


0.001 - • 


FSRQ 
BL Lac 
3C 273 
Mkn 421 
HSP BL Lacs 
Unidentified 


• : • 

* • • 


10 


36 


•. V* *• 


1 o 37 10 38 10 

Jet power (W) 

(a) 


39 


10 


40 


o 

Cl 


>< 100 - 


o 

C3) 

CD 

i. 

O) 

g 

H—' 

05 

i. 

Q 

CD 

O 

O 

05 

H— 

o 

SI 

- 4 —* 

C3) 

c 

CD 


10 - 


1 - 


0.1 - 


FSRQ 
BL Lac 
3C 273 
Mkn 421 
HSP BL Lacs 
Unidentified 


x 

x 


x 

x 


• • 


•• 


10 


20 


40 


Maximum bulk Lorentz factor 

(b) 

Figure 3. The results of fitting our model to all 38 blazars with simultaneous multiwavelength spectra and redshifts from Abdo et al. 2010 and the 4 quasi- 
simultaneous HSP BL Lac SEDs from Padovani et al. 2012. Figure a shows the radius of the jet transition region plotted against the jet power for all fits 
(the transition region is where the jet transitions from parabolic to conical, stops accelerating and first approaches equipartition due to strong non-thermal 
particle acceleration). We see a clear correlation between the jet power and radius of the transition region as suggested in our previous work. BL Lacs have 
systematically smaller radius transition regions than FSRQs, with an approximately linear correlation between the two as illustrated by the dashed lines. We 
have highlighted the two outliers to the relation, Markarian 421 and 3C273. We argue that these two blazars are atypical since Mkn421 is one of the closest 
known blazars and 3C273 is the one of the most powerful (and seems to be misaligned). Figure b shows the maximum bulk Lorentz factor plotted against the 
distance of the transition region from the black hole (the ratio of transition region radius to length is fixed at ~ 1 : 80). We see that the bulk Lorentz factor is 
systematically larger for FSRQs than BL Lacs and increases with the distance over which the jet accelerates as expected for steady magnetic acceleration. The 
relationship between the length of the accelerating region and the maximum bulk Lorentz factor is approximately 7 max oc cc^/ 4 , or equivalently in terms of 
the transition region radius 7 max oc 
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3 RESULTS 

Our model is able to fit to all 38 Fermi blazars with simultane¬ 
ous multiwavelength spectra and redshifts from Abdo et al. (2010) 
and the 4 high power, high synchrotron peak frequency (HSP) BL 
Lacs from Padovani et al. (2012) with unprecedented precision. The 
spectral fits are shown in Figure 6 and their parameters in Figure 7. 
One of the main achievements of this work is to show that an ex¬ 
tended fluid model can fit to both the radio and high energy emis¬ 
sion simultaneously. Fitting to the radio data simultaneously with 
all other wavelengths allows tight constraints to be placed upon the 
radius of the brightest synchrotron emitting region (in our model 
this is the transition region where the jet first becomes conical, stops 
accelerating and comes into equipartition. see Fig. 1). We have sum¬ 
marised our main results in Figure 5. 


3.1 Constraining the radius of the transition region 

The radius of the brightest emitting synchrotron region is tightly 
constrained by modelling the break frequency where the syn¬ 
chrotron emission goes from being self-absorbed and close to flat 
in flux at low frequencies, to optically thin at high frequencies (the 
synchrotron break has also been used to constrain the properties of 
a small number of X-ray binary systems with observed jet spectra 
e.g. Russell et al. 2013). This break at low frequencies ~ 10~ 4 eV 
(~ 20GHz) is clearly seen in many of the spectra e.g. J0423, J0531, 
J1057, J1159, J1221, J1256 and J1512 in Fig. 6 (for a simplified 
derivation of the break frequency and an explicit plot illustrating 
the sensitive dependence on the jet radius see equation 8 and Fig. 
2 in Potter & Cotter 2013b) . The determined radius of the tran¬ 
sition region is shown for all 42 blazars in Figure 3a. There is a 
clear linear correlation between the radius of the transition region 
and the jet power. This is significant because it shows that there 
is a simple underlying relation between the structural properties of 
jets and their power. In this work we have assumed that the transi¬ 
tion region is the location where the jet stops accelerating and first 
transitions front a parabolic to a conical jet based on the results of 
GRMHD simulations. This means that the initial power of the jet 
determines the radius at which it will stop accelerating and start 
to brightly emit synchrotron radiation. This result can be used to 
inform and constrain the acceleration mechanism of the jet. Cur¬ 
rently, simulations tend to find that jets stop accelerating at roughly 
10‘Vs (McKinney 2006 and Komissarov et al. 2007), however, the 
transition from parabolic to conical in M87 occurs at 10 5 r s . These 
observations combined with our findings suggest that real jets can 
accelerate over much larger distances than currently found in sim¬ 
ulations. This difference could be a result of the interaction of the 
jet with its host environment and its initial mass-loading combined 
with the effects of radiative energy losses and drag, which are diffi¬ 
cult to include properly in simulations. 


3.2 A relation between the length of the accelerating region 
and jet speed 

Figure 3b shows the length of the accelerating region of the jet (our 
jet shape is linearly scaled from M87 and so the radius of the transi¬ 
tion region is « 1 /80 of the jet length at that point) plotted against 
the maximum bulk Lorentz factor at this radius. There is a clear cor¬ 
relation between the length of the accelerating region and the max¬ 
imum jet speed, with evidence for bimodality. This suggests that a 
steady acceleration process is operating in the jet. The longer the 
distance over which the jet accelerates, the higher its bulk Lorentz 


factor becomes. This is consistent with steady magnetic accelera¬ 
tion in the parabolic section of the jet, however, it also shows that 
the final Lorentz factor of the jet depends on its initial jet power. 
We also see a clear association of FSRQs with fast jets and BL Lacs 
with slower jets. In simulations the maximum bulk Lorentz factor 
depends primarily on the initial mass loading or equipartition ratio 
of the magnetic field to particle energies. This is because the jets 
only accelerate efficiently up to equipartition, where approximately 
half of the initial magnetic energy is converted into bulk kinetic en¬ 
ergy in the particles i.e. 


(ML 


E p 

~Eb' 


7max 


1 + (ML 

2 (ML 


(25) 


where (ml is the mass loading fraction and 7 max is the terminal 
maximum bulk Lorentz factor. From Fig. 3b we find a relationship 
7max oc R-f 4 between the radius of the transition region and the 
maximum bulk Lorentz factor of the jet. This is a weaker depen¬ 
dence of the maximum bulk Lorentz factor with jet radius than 
typically found in simulations 7max oc R 1 ' 2 McKinney (2006); 
Komissarov et al. (2007), although often the chosen prescription 
for the confining wind or rigid jet boundary will determine this re¬ 
lationship in simulations. This difference could be due to the en¬ 
trainment of material surrounding the jet or radiative drag on the 
jet plasma as electrons lose bulk momentum by inverse-Compton 
scattering external radiation fields. 

This is the first attempt to constrain the relation between the 
length of the accelerating region of the jet and the maximum bulk 
Lorentz factor using a large sample of blazars, so these results are 
very useful to compare the results of simulations and observations. 
One potential problem in interpreting our results is our assumption 
that the brightest synchrotron emitting region occurs when the jet 
first comes into equipartition and stops accelerating. This scenario 
seems to us the most likely, especially given the close proximity of 
the bright feature HST1 in M87 with the transition from parabolic 
to conical Asada & Nakamura (2012). However, it is possible that 
the region where substantial particle acceleration occurs is not close 
to the transition between parabolic and conical. It is worth empha¬ 
sising that the tight constraints our model places on the radius of 
the brightest synchrotron emitting region are independent of the as¬ 
sumptions we have made about the accelerating region of the jet. 
Irrespective of the validity of our individual assumptions our results 
show the importance of using extended fluid jet emission models to 
understand the relationships between basic structural and dynami¬ 
cal quantities, since these are crucial to informing, constraining and 
testing our theoretical models and the results of GRMHD simula¬ 
tions. 


3.3 Evidence for a bimodal accretion rate and AGN 
unification 

In order to understand and interpret our results in an astrophysi- 
cal context it is necessary to make assumptions about the black 
hole mass distribution in the blazar population. Although black hole 
masses are notoriously difficult to determine accurately, evidence 
suggests that the black hole masses of BL Lacs and FSRQs (gener¬ 
ally associated with low and high power jets) are similar, with the 
typical variation in the black hole masses being roughly an order of 
magnitude (Shaw et al. 2013). These results suggest that there is no 
large systematic variation in the black hole masses of blazars and 
therefore that the black hole mass is unlikely to be driving the corre¬ 
lations shown in Fig. 3, which exist over many orders of magnitude. 
In order to interpret our results in terms of the Eddington luminos- 
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(a) (b) 

Figure 4. The distribution of fractional Eddington luminosity and the distance to the jet transition region under the assumption of a constant black hole mass of 
5 X 10* ATq for all 37 blazars with simultaneous multiwavelength Fermi observations and classifications. Figure a shows a clear dichotomy in the fractional 
Eddington luminosity between BL Lacs and FSRQs, suggestive of two accretion modes. Figure b shows that the length of the accelerating parabolic region of 
the jet, xt, extends out to larger distances in FSRQs than BL Lacs. This could be responsible for the different maximum bulk Lorentz factors shown in Figure 
3b. 


ity we shall therefore assume a fiducial blazar black hole mass of 
5 x fO s M 0 (Shaw et al. 2012, 2013). Under this assumption we 
find evidence for a bimodal distribution of accretion rates in the 
blazar population (Fig. 4). We find that FSRQs are high power jets 
close to the Eddington luminosity whilst BL Lacs are much lower 
power jets, consistent with previous work Ghisellini et al. (2011). 
These results agree remarkably well with the results of AGN lu¬ 
minosity distributions (see Fig. 14 in Mingo et al. 2014) with FS¬ 
RQs corresponding to high excitation radio galaxies (HERGs) with 
0.05Ledd — 1-Ledd an d BL Lacs corresponding to low excita¬ 
tion radio galaxies (LERGs) with 2 x 10~ 4 Ledd — 0.05-Ledd- 
This close agreement between our model and the observations of 
AGN acts as an important, independent verification of the results 
and assumptions of our model. This strongly suggests that there is 
a correspondence between FSRQs and HERGs (closely associated 
with FRII sources) and BL Lacs and LERGs (closely associated 
with FRI sources) supporting the idea of AGN unification. 

Under the same assumption of a fiducial black hole mass we 
find a bimodal distribution of the length of the accelerating region 
of a jet and its power in terms of the Eddington luminosity. This 
suggests that for a given black hole mass more powerful jets will 
accelerate over larger distances and reach higher terminal veloci¬ 
ties. This could be interpreted in terms of a more powerful jet trav¬ 
elling a larger distance before being disrupted by its environment 
(see e.g. Bicknell 1995). We also find a close association of FSRQs 
with fast jets and BL Lacs with slower jets. Together Figures 4a and 
4b provide strong evidence for the existence of bimodal accretion 
rates and association of FSRQs with powerful, fast FRII type jets 
and BL Lacs with less powerful and slower FRI type jets. 

3.4 A physical understanding of the blazar sequence 

The results of this work allow for a physical interpretation of the 
blazar sequence. The blazar sequence, originally discovered by Fos- 
sati et al. (1998), describes the combination of an observed anti¬ 
correlation between the jet power and the synchrotron peak fre¬ 
quency in the blazar population, and a correlation between the jet 
power and the Compton dominance (the ratio of synchrotron to 
inverse-Compton luminosity). This result has been the subject of 


considerable discussion and speculation in recent years (e.g. Ghis¬ 
ellini et al. 1998; Ghisellini & Tavecchio 2008; Maraschi et al. 
2008; Meyer et al. 2011), however, our work is the first systematic 
attempt to simultaneously model both the radio and high energy 
emission of a large sample of blazars to understand their structural 
and dynamical properties and their relation to the blazar sequence. 
We find that BL Lacs are low power blazars, with jets that accelerate 
over shorter distances and reach lower maximum Lorentz factors 
than more powerful FSRQs. This has several implications for their 
spectra: firstly the brightest synchrotron emitting region in BL Lacs 
is smaller in radius than in FSRQs and this means that they have 
larger magnetic field strengths. This is easily calculated by equat¬ 
ing the rest frame magnetic energy contained in a thin slab of the 
jet plasma in equipartition, to half the total energy of the slab 


b' t ocASL 

3 2 no 2 Tmaxfir 


(26) 


where a subscript T indicates the value of a quantity at the transi¬ 
tion region, Wj is the lab frame jet power and l s is one light-second 
l s = cx ls. Since BL Lacs have both lower maximum bulk Lorentz 
factors and smaller transition region radii than FSRQs (from Fig. 
3a we find Wj oc Rt) the magnetic field strength at the transi¬ 
tion region is greater (whilst the jet power is smaller, this factor is 
suppressed by a square-root). The synchrotron peak frequency is 
directly proportional to the magnetic field strength and so we find 
that as the jet power and jet radius increase the synchrotron peak 
frequency decreases. This offers a physical explanation for the an¬ 
ticorrelation between jet power and peak synchrotron frequency. 

The inverse-Compton emission from BL Lacs is dominated by 
synchrotron self-Compton (SSC) emission due to density of syn¬ 
chrotron photons in the compact, bright transition region with a 
large magnetic field strength. We find that the power emitted by 
SSC is naturally limited to being comparable to, or smaller than, 
the power emitted by synchrotron emission and this prevents SSC 
dominated sources (i.e. most BL Lacs) from becoming Compton- 
dominant. This is because the highest energy electrons are both 
responsible for emitting the majority of the synchrotron power 
(oc Es) and SSC power (oc EsU sync h ~ E^). This means that as 
the spectrum becomes Compton-dominant, the SSC radiative cool- 
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BL Lacs - Small radius transition region 
with large magnetic field strength 
leading to high peak synchrotron 
frequency. Gamma rays dominated by 
SSC emitted within transition region. 


BL Lacs: y bu[k ~ 6-20 

0.1/10pc- 5 x10 38 W (~0.05I_edd) 

FSRQs: Xftui/c ~20 

2/200pc 10 4 °W (~1 Ledd) 


R t oc Jet Power 

Max. bulk Lorentz factor oc R T 1/4 



Figure 5. A summary of our main results. The reader may find it useful to refer to the jet schematic in Fig. 1 to understand the jet geometry and basic 
assumptions of our model. 


ing lifetimes for these high energy electrons decreases rapidly and 
so it is difficult to replenish this high energy population and main¬ 
tain Compton-dominance via SSC (unless the electron spectrum, 
N e (E e ) oc Ey a , has a very low spectral index, i.e. a < 1, so that 
the majority of the energy in the electron population is concentrated 
in the highest energy electrons). 

In FSRQs we find that the radius of the brightest synchrotron 
emitting region, constrained by fitting the synchrotron break, is too 
large, and the magnetic field strength too small, for substantial SSC 
emission to occur. This is because SSC emission requires a large 
energy density of synchrotron seed photons to Compton up-scatter 
and the power emitted by synchrotron radiation is proportional the 
the square of the magnetic field strength. For FSRQs we find the 
radius of the transition region to be between 0.1pc-2pc and the 
corresponding distance to be « 10pc-200pc. At these large dis¬ 
tances which are well outside of the BLR (7?blr < lpc, e.g. Kaspi 
et al. 2005) and outside the dusty toms (77dt < lOpc, e.g. Elitzur 
2006 and Nenkova et al. 2008), we find the high energy emission 
is well fitted by scattering CMB photons. Since inverse-Compton 
scattering of external photon sources is generally required to fit to 
most low-peak frequency, Compton-dominant FSRQ spectra, it is 
no surprise that most previous studies have required the brightest 
emitting region to occur within the BLR or dusty toms, since these 
were the only external photon sources which were modelled (e.g. 
Blazejowski et al. (2000), Ghisellini & Tavecchio (2008), Sikora 
et al. (2009). Tavecchio et al. (2011) and Joshi et al. (2014)). How¬ 
ever, crucially, these works neglected to model the radio observa¬ 
tions. We find that by simultaneously modelling both the radio and 
high energy emission we are able to constrain the brightest region 
in FSRQs, to occur at much larger distances than the BLR or dusty 
torus. 

3.5 Limitations of a time-independent model 

In this and our preceding papers we have developed the first rela¬ 
tivistic fluid jet model which calculates the emission from all sec¬ 


tions of the jet. We decided to initially study a time-independent 
model in order to understand the average properties of the jet fluid 
flow before adding time-dependent perturbations. Whilst this has 
the benefit of allowing us to understand correlations between the 
long-term average jet properties it does not allow us to investigate 
the dramatic short-term variability of blazars. Bearing in mind the 
very short timescales for gamma-ray flares (e.g. Albert et al. 2007, 
Aharonian et al. 2007 and Abdo et al. 2010) and the range in flaring 
behaviour, we expect that whilst the long-term steady gamma-ray 
emission seems to be well fitted by scattering CMB photons at large 
distances, it seems very unlikely that flares at such large distances 
could be responsible for the very short timescale gamma-ray vari¬ 
ability observed due to the long radiative cooling and light-crossing 
timescales. We therefore expect that short timescale flaring events 
are likely to occur at much shorter distances along the jet, probably 
within the parabolic accelerating region. We shall investigate vari¬ 
ability and flaring by adding time-dependent perturbations to the 
steady fluid flow in the near future. 


4 CONCLUSION 

In this work we have used one of the most realistic jet emission 
models to date to fit with unprecedented accuracy to the spectra of 
the entire sample of Fermi blazars with simultaneous multiwave¬ 
length observations and redshifts from Abdo et al. (2010). From 
this we have been able to place constraints on the basic structure 
and dynamics of blazar jets. Our main results are as follows: 

• We find a linear correlation between the radius of the brightest 
synchrotron emission region and the jet power, which we constrain 
via the optically thick to thin synchrotron break. In most FSRQs 
this radius (>0.1pc) is large enough to place this region, containing 
the most intense particle acceleration, well beyond the broad line 
region and dusty torus. 

• We find a correlation between the final radius of the accelerat- 
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ing region of the jet and the maximum bulk Lorentz factor of the jet 
given approximately by 7buik oc T?^/ 4 . 

• There is a bimodal distribution of jet power in the blazar pop¬ 
ulation with FSRQs having large jet powers (5 x 10 38 W — 10 40 W) 
and BL Lacs (5 x 10 36 W — 5 x 10 38 W). Assuming a fiducial black 
hole mass of 5 x 10 8 Mq we find excellent agreement with the 
distribution of AGN luminosities from Mingo et al. (2014), with a 
close correspondence of FSRQs with high excitation radio galaxies 
(FIERGs, 0.05Ledd — ILedd) and BL Lacs with Low excitation 
radio galaxies (LERGs, 2 x 10 _4 Ledd — 0.05Ledd) • 

• There is a bimodal distribution of maximum bulk Lorentz fac¬ 
tors in the blazar population with FSRQs having large bulk Lorentz 
factors (20-50) and BL Lacs lower bulk Lorentz factors (6-20). To¬ 
gether with the differences in jet power, this provides evidence for 
AGN unification with FSRQs corresponding to powerful, fast FRII 
type jets and BL Lacs corresponding to weaker, slower FRI type 
jets. 

• FSRQs have larger maximum bulk Lorentz factors and a larger 
radius when the jet first comes into equipartition and intense parti¬ 
cle acceleration occurs. This means they have lower magnetic field 
strengths in this bright region, which leads to sub-dominant SSC 
emission and low peak synchrotron emission. We find that for FS¬ 
RQs, their large bulk Lorentz factors, maintained to large distances, 
allow the Compton-dominance and low inverse-Compton peak fre¬ 
quency to be explained via inverse-Compton scattering of CMB 
photons. 

• BL Lacs have lower bulk Lorentz factors and smaller radius 
jets at the location of the brightest synchrotron emission, leading 
to higher magnetic field strengths in this region than FSRQs. This 
means that the peak synchrotron emission from BL Lacs occurs at 
higher frequencies and SSC dominates the inverse-Compton emis¬ 
sion. These dynamical and structural constraints allow an under¬ 
standing of the physics underlying the famous blazar sequence. 

These results are summarised in Fig. 5 and demonstrate the 
merit of using realistic, extended fluid jet emission models to fit to 
blazar spectra and constrain their basic properties. 
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Figure 6 . The fits to the SEDs of all 38 Fermi blazars with simultaneous multiwavelength observations and redshifts from Abdo et al. 2010, and the 4 non- 
simultaneous spectra of the high synchrotron peak frequency FSRQs from Padovani et al. 2010 (J0035, J0537, J0630 and J1312). The model is able to fit very 
well to all of the spectra including radio observations. Simultaneous (and quasi-simultaneous) data points are shown as solid circles whilst archival data are 
included as crosses. Due to computational time constraints fits were performed by eye. 
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Figure 7. Fitting parameters for all 42 blazars used in Figs. 6 and 3. Columns from left to right: Fermi blazar name (RA only), initial jet power, jet length, 
electron cutoff energy, injected electron power law index, jet misalignment angle, bulk Lorentz factor at transition region, bulk Lorentz factor at end of jet, 
accretion disc luminosity, radius of transition region, magnetic field strength at transition region and effective black hole mass (if the jet transition region were 
to occur at 10 5 r s as in M87) also used as the black hole mass for the thin accretion disc fit. 
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